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ABSTRACT 

Context. Surveys of the binary populations in the solar neighbourhood have shown that the periods of G- and M-type stars are log- 
normally distributed in the range from 0.1 - 10" days. However, observations of young binary populations in various star forming 
regions suggest a log-uniform distribution. Clearly some process(es) must be responsible for this change of the period distribution 
over time. Most stars form in star clusters, so it is here that the(se) process(es) take place. 

Aims. In dense young clusters two important dynamical processes occur: i) the gas-induced orbital decay of embedded binary systems 
and ii) the destruction of soft binaries in three-body interactions. The emphasis in this work is on orbital decay as its influence on the 
binary distribution in clustered environments has been largely neglected so far. 

/Wefhods. We performed Monte-Carlo simulations of binary populations to model the process of orbital decay due to friction with 
the gas. In addition, the destruction of soft binaries in young dense star clusters was simulated using nbody modelling of binary 
populations. 

Results. It is known that the cluster dynamics destroys the number of wide binaries, but leaves short-period binaries basically undis- 
turbed. Here it is demonstrated that this result is also valid for a initially log-uniform period binary distribution. By contrast the 
process of orbital decay significantly reduces the number and changes the properties of short-period binaries, but leaves wide binaries 
largely uneffected. Until now it was unclear whether the short period distiribution of the field is unaltered since its formation. It is 
shown here, that if any alteration took place, then orbital decay is a prime candidate for this task. In combination the dynamics of 
these two processes, convert evan an initial log-uniform distribution to a log-normal period distribution. The probability is 94% that 
the evolved period distribution and the observed period distribution have been sampled from the same parent distribution. 
Conclusions. Our results provide a new picture for the development of the field binary population: Binaries can be formed as a result 
of the star-formation process in star clusters with periods that are sampled from the log-uniform distribution. As the cluster evolves, 
short-period binaries are merged to single stars by the gas-induced orbital decay while the dynamical evolution in the cluster destroys 
wide binaries. The combination of these two equally important processes reshapes a initial log-uniform period distribution to the 
log-normal period distribution, that is observed in the field. 

Key words. Keywords should be given 



1. Introduction 

The most important observations that shaped our current pic- 
ture of the binary field population were already performed in the 
early 1990's. Duquennoy & Mayor ( 1991 ) determined the binary 
frequency of G-type stars in the solar neighbourhood to be about 
61% and the period distribution to follow a log-normal distribu- 
tion. 



/(logP) = Cexp 



-(log P- log P)^ 



(1) 



over the period range ~ 10 ' - 10"d where P is the period in 
days (log P - 4.8 = 1 72yr, o-jo^p - 2.3) and C a normalisation 



constant. A year later ' Fischer & Marcy| ( |1992| l analysed the bi- 
nary properties of M-dwarfs in the separation range 0.04- 10"* AU 
corresponding to periods in the range ~ 10 2 - lO^'yr and found 
that the period distribution of M-dwarfs in the solar neighbour- 
hood is also log-normally distributed with a peak between 9yr 
and 270yr - nearly identical to the findings for the G-dwarfs. 
Although the observed binary frequency of M-dwarfs with about 
33% is lower than that of G-dwarfs, the properties of the G and 
M-dwarf binary populations seem very similar. More recent ob- 



servations by Raghavan et al. ( 2010| l confirm this picture. They 
determined the binary properties of nearby (d < 25pc) solar- 
type stars (~ F6 - K3) and find the fraction of binary stars to 
be 58% with a log-normal period distribution (logP - 5.03 and 
o" logf - 2.28). 

With an average age of a few Gyrs the field population con- 
stitutes an dynamically evolved state. Therefore, dynamical pro- 
cesses have most likely changed the binary properties since the 
formation of these stars and the current properties differ signif- 
icantly from the primordial state. In order to understand the bi- 
nary formation process it is not sufficient to know the properties 
of the (old) field population but those of the primordial binary 
population. 



Connelley et al. ( 2008 1 determined the binary properties of 
the very young populations in Taurus, Ophiuchus and the Orion 
star-forming regions (excluding the much denser Orion Nebula 
Cluster). They found that their observed peiiod distribution can 
be fitted by a log-uniform distribution, differing significantly 
from the log-normal distribution in the field (see also |Kraus &| 
Hillenbrand ( 2007 )). In these sparse young star forming regions 
it is improbable that dynamical evolution altered the period dis- 
tributions in the short period since the stars were formed (see 
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e.g. |Kroupa & Bouvier||2003| l. So it can be assumed that their 
properties match the initial conditions in general. 

In denser regions indications for an originally log-uniform 
distribution are found, too. For example, HST observations by 



Reipurth et al. ( 2007 1 of binaries in the Orion Nebula Cluster 
(ONC) demonstrate, that the semi-major axis distribution de- 
viates from the log-normal distribution and is closer to a log- 
uniform distribution. So it seems that older binary populations 
have a log-normal period distribution while the primordial dis- 
tribution is likely to be log-uniform. 

From the theoretical side the binary frequency is the best 
studied binary property (e.g. Kroupa & Burkert 2001; Kroupa| 
& Bouvier.2003, . Ki'oupa et al..2001, , Pfalzner & Olczak,2007f 



Marks et al.|201 l[|Marks & Kroupa|20TT| i7lt seems that also the 
evolution of the binary frequency in the ONC depends on the 
initial binary frequency, the evolution of the binary population 
does not ( [Kaczmarek et al.|201 1 ). However, the binary popula- 
tion is not only described by the binary frequency but also by the 
period, mass-ration and eccentricity distribution. 

Starting with the work by Heggie ( 1975 ), it has been realised 
that binaries become dynamically destructed via three- and four- 
body interactions. Generally wide binaries are more affected by 
dynamical destruction than close ones. The existence of bina- 
ries wider than 10^ AU in the field still poses an open question 
regarding their origin ( Parker et al.|2009|. 

Performing N-body simulations [Kroupa (T995a]b| showed, 
that to reproduce the observed log-normal distribution of the 
field, the initial number of wide binaries has to be significantly 
higher than observed, if all binaries are exposed to dynamical 
evolution in a star cluster This rising distribution was obtained 
by inverse dynamical population synthesis, inverting the effects 
of dynamical destruction on the period distribution (see dashed 
line inFigJTjl. 

Parker et al. ( 2009 | l finds that binaries with periods exceeding 



lO"' AU cannot survive in these clusters he investigated. However, 
binaries with semi-major axes exceeding lO^'AU are observed in 
the field. [Kouwenhoven et al. ( 2010 1 suggest, that these observed 
wide binaries could form during the cluster dissolution. Another 
approach would be a different cluster type. 

Here we investigate how the initial period distribution would 
evolve in its natal environment considering the two different dy- 
namical processes that inevitably affect binaries in a cluster: i) 
gas-induced orbital decay and ii) dynamical destruction of bina- 
ries in encounters. 

So far the influence of orbital decay has been only investi- 
gated for isolated binary systems Stabler ( 2010| l. Orbital decay 
takes place in the earliest stages of star formation, when the stars 
are still embedded in the gas. Here the interaction of binaries 
with the gas leads to the excitation of waves in the surround- 
ing gas. The energy transfer from the binary to the gas leads to 
orbital decay ( [Stabler [2010] ). Here we investigated how orbital 
decay changes the period distribution in a cluster environment. 

This paper is structured in the following way: in Section|2]the 
cluster setup is explained. Section |3] describes the orbital decay 
process in the gas embedded phase. We demonstrate in Section]?] 
that a combination of the orbital decay with the dynamical de- 
struction ( [Kaczmarek et al. ] [2011| l naturally transforms a log- 
normal birth period distribution into a log-uniform period dis- 
tribution. Our conclusions mark the last section. 



2. Cluster Setup 

The processes of orbital decay and dynamical destruction take 
inevitably place as soon as the binary population forms in a 
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Fig. 1. Schematic picture of possible initial period distributions. 
The thick grey curve displays the log-normal fit to the data of 
[Raghavan et al.[ ( [2010| l, the thick dashed line the |Kroupa] ( [1995a| l 
and the solid line a log-uniform period distribution, that is ob- 
served in many young star forming regions (see text for refer- 
ences). 



dense star cluster environment. Both result from interactions be- 
tween stars in a dense environment. As most stars form in cluster 
environments ( Lada & Lada 2003 ), it can be anticipated that the 
field binary population we observe today was - at least to same 
extend - shaped by these processes. 

In this study we chose the Orion Nebula Cluster (ONC) as 
model cluster It is probably a typical star forming region and 
has the advantage that many parameters both of the cluster struc- 
ture and the binary population are well observed (e.g. [Jones & 



Walker'1988';'Prosseretal.'1994';'Hinenbrand'1997rHilienbrand 
& Hartmann 1998; O'Dell et al. 2009; Preibisch et al. 1999] 
Kohler et al. 2006 ; Rei purth et al.f2007] ). The ra dial density pro- 
file declines in the outer parts as Ppresent '"^ (Jjones & Walker 



1988 



Hillenbrand 



1997[) and is much flatter (pcore °^ r ^ within 



^core ~ 0.2pc. The dynamical evolution of the cluster signifi- 
cantly changes the stellar density profile during the first Myr, 
the estimated age of the ONC (Hillenbrand 1997| l. Hence we 
adopted the following initial density profile 



[ 0, 



core 



Q<r<R, 
RcoK <r<R 
R < r < oo. 



(2) 



Here we choose Rcors = 0.2 pc and a cluster size of R = 2.5 pc 
to model the ONC. The used model starts from a situation were 
all stars are already formed. In dynamical equilibrium the initial 
radial density profile (Eq. ]2|l evolves towards the currently ob- 
served one over a time span of about IMyr ( [Olczak et al.|201()) l. 
However, other initial conditions could result in the same profile 
( [Kroupa|2000]|Amson et al.|2009| l. 

To mimic the often observed mass-segregation in young 
dense clusters, the most massive binary system is initially placed 
at the cluster centre and the three next most massive stars at ran- 
dom positions within a sphere of radius R = 0.6 Rhm around 
the cluster centre, where Rhm is the cluster half-mass radius. 

T99^ 



In this we follow the suggestions by Bonnell & Davies 



who found that the observed mass segregation (Hillenbrand & 
[Hartmann" 1998) of the ONC is unlikely to be the result of the 
dynamical evolution of the cluster but has to be primordial, if 
the cluster is initially in virial equilibrium. Note, that if the clus- 
ter would be in a subvirial state, probably no primordial mass 
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segregation would be required ( [Allison et al.|2009t|Olczak et al.| 
[2UTTI 1. 

All stellar systems are initially created as binaries with sys- 
tem masses sampled from the initial mass function (IMF) sug- 
gested by Kroupa ( 2007^ Eqs. 19 & 20). We limited our system 



mass Msys = M^i + Ms2 to [O.OSMq, 50Mq], where M.51 and 
Ms2 are the masses of the stars constituting a binary system. The 
upper mass limit corresponds to the mass of the most massive bi- 
nary system of the ONC. The lower limit is the hydrogen burning 
limit, and thus the margin between sub-stellar and stellar objects. 

Ideally we would like to use observed initial properties of 
binaries in very young embedded clusters. We will see that em- 
bedded clusters younger than 100,000 yr would be ideal, because 
at such a young age the binary population would be close to its 
primordial state. However, it is very difficult to observe clusters 
at such a young age. Observations of binaries in dense embed- 
ded clusters face a number of observational difficulties like ex- 
tinction, crowding, resolution etc.. In addition, the number of 
stars would be very low, making statistical statements impossi- 
ble. As their is an age spread in forming clusters (1-3 Myr) it 
would be a possibility to use just the youngest stars in a form- 
ing cluster to deduce the primordial state. However, age deter- 
mination for such young stars is problematic. Alternatively, one 
can assume that low density clusters are relatively unaffected by 
these processes described here and are therefore close to their 
primordial state. Observations of the primordial state not be- 
ing available, we use the work by Kouwenhoven et al. ( [2007| l 
based on observations of the low-density Scorpius OBII cluster. 
Choosing the semi-major axis and mass ratio distribution from 



an OB association with an age of (5-20) Myrs (Kouwenhoven 
[et al. 2007 ), as Scorpius OBII, seems at first glance unusual, as 
at that age the binary population might have already been pro- 
cessed. However, the low stellar density of Scorpius OBII com- 
bined with the lack of surrounding gas speaks for the cluster 
having neither been processed by orbital decay nor by dynami- 
cal destruction, as mentioned before. 

The primordial semi-major axis distribution was chosen 
to be log-uniform /«(«) oc dN/da cx a"' with a € 
[0.02 AU, 10000 AU]. As primordial mass ratio distribution we 
started with fg(q) oc ^"-^ and q - MsilMsi € [0, 1]. 

Another justification for our choice of the initial semi-major 
axis distribution, is that the log-uniform distribution is a simple 
straight forward distribution, that requires strong processing in 
order to obtain the log-normal field distribution. This distribution 
was used here to test, if even this extreme assumption leads to the 
log-normal distribution observed in the field. This does not rule 
out, that one obtains a similar result as described in the following 
with different initial conditions, containing fewer wide and/or 
close binaries. 

In the following the above described model will be used to 
study the effect of the gas-induced orbital decay of binaries in 
star clusters. 

3. Gas-induced Orbital Decay 

In the early embedded phases clusters consist of the already 
formed stars and a large gas component from which further stars 
potentially form. These stars keep being embedded in their natal 
gas cloud until it is removed by strong winds, radiation of the 
massive stars and supernova explosions. During this embedded 
phase, stars and binary systems experience dynamical friction 
with the ambient gas. Binary systems induce spiral waves in the 
gas leading to energy and angular momentum loss resulting in 
a shrinking of the binary orbit and possibly merging of the two 



stars . jStahlerj ( [20 1 0) l derived an analytic expression for the tem- 
poral development of the separation a,„, for a single isolated bi- 
nary system on a circular orbit. The separation diminishes due 
to gas-induced orbital decay with time as 



flO 1 



(3) 



where aq is the initial binary separation and tc the so-called coa 
lescence time, which is given by 

15 d+qfc] 



tr = 



flo 



32;r q poG^M%^' 



(4) 



where G is the gravitational constant, c, the sound speed, q the 
mass ratio and Msys the system mass of the binary system. 

Stahler (2010) assumed in his derivation that gas interactions 
close to the stars themselves can be neglected. The binary sys- 
tem, represented by an oscillating gravitational potential, torques 
the nearby gas and produces outgoing acoustic waves. These 
waves transport angular momentum from the binary to the sur- 
rounding gas. Thus the orbit of the binary decays as long as the 
gas density is high enough. Eq.[3]is only valid for radial distances 



< r 



GMsys""" 

2 c? 



Ill AU 



rMsys™"^ 



2 km s 



(5) 

between the two stars constituting the binary system. 

To analyse the effect of gas-induced orbital decay of binary 
systems on the period distribution, we applied Eq.|3]to the entire 
population of a young star cluster, for the example of the ONC. 
Starting from a cluster with initial conditions as given in Sec.|2 
the temporal development of the binary population is followec 
for different gas density models. 

Only if a binary system fulfills the restriction given by Eq. |5] 
the orbital decay is calculated, otherwise the orbit remains unal- 
tered. This means, for example, for a system mass of 0.08 M0, 
1 Mq and 50 M© the orbital decay will only be calculated for bi- 
nary systems with a separation r < 8.88 AU, r < 111 AU and 
r < 5550 AU. Like in Stahler's work modelling the orbital de- 
cay the simplification of treating all orbits as circular has been 
adopted. Future work should include elliptical orbits as well. 

For the gas-induced decay, the gas density distribution in the 
cluster is of vital importance. Here we assume an isothermal gas 
density which follows the stellar density distribution. To prevent 
the distribution from diverging at the centre, the density is kept 
constant at a value p^ax inside the cluster core radius Rcok = 
0.2pc. Outside this area the density decreases isothermally. Thus 
the isothermal gas density distribution can be described by the 
following equation 



Pgas('") = Pn 



1, 



(^core/^) 5 ^core f ^cluster 



(6) 



For an isothermal density distribution Pgas('") = c JlnGr 



(e.g. Binney & Tremaine 1987 Eq. 4.123) the sound speed is 
given by 



(7) 



Although the density distribution (Eq. |6]l is not isothermal 
over the whole parameter range, we approximate the sound 
speed as 



^J2nGp, 



max ' core • 



(8) 
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The resulting sound speeds for distributions with p,„„v e 
[10"^, lQ^]cm^^ lie between the sound speed of an ideal gas at 
lOK, Cj - Q.2km/s, and the observed value in infrared dark 
clouds ( |Sridharan eTaL|[2005] ) q = [1.0 - 2.2]km/s. A even 
better approach would be a time dependent density p{r, t) and 
sound speed Cs{r, t), which hasn't been observed yet. Deviations 
of the sound speed manifest themselves in the coalescence time 
(Eq.g. 

If the separation becomes smaller than the radius of both 
stars fltot < Rsi + (assuming a main-sequence mass-radius 
relationship (|Binney & Merrifield 1998 p. 110)), 



R 



.sl/2 



M,i/2 



0.7 



a binary system is treated as a 'merged' system. This is clearly 
a strong simplification, because the merging of stars is a much 
more complex process that depends on a variety of parameters 
such as the presence of a surrounding disc, the eccentricity of the 
orbit, mass transfer and the conditions of the molecular cloud in 
which the stars are embedded. However, including all these pro- 
cesses goes beyond the aim of the current paper. Here, all pro- 
cesses at these small distances are excluded under the assump- 
tion that these two stars might possibly merge at such small dis- 
tances. 

First it was analysed how the binary population parameters 
change due to orbital decay in general. Therefore, a series of 
simulations with three different gas densities in the range of 
[lO'^cm"-', lO^cm"-'] were performed corresponding to moder- 
ate gas densities observed in star forming regions (Padmanabhan 
[200 1 1 . For the three densities lO'^cm ^ 10^ cm^^ and 10'' cm"^' 
the sound speeds are 0.15 km/s, 0.49 km/s and 1.54 km/s, re- 
spectively (Eq.[8]). The results summarised in Table [T] show the 
merger rates for all combinations of gas densities and sound 
speeds. As to be expected, the percentage of possibly merged 
binary systems after 1 Myr (Table [TJ depends on the maximum 
density of the isothermal density distribution (Eq.[6]l. The value 
of the overall reduction of binaries could vary for a constant 
sound speed by a factor of two to four if the density would 
change by one order of magnitude. Keeping the density constant, 
one can see, that the percentage of merged systems is sensitive 
to the choice of sound speed. 

Additionally, the degree of orbital decay depends on the gas 
density distribution. Here only the case of a temporary con- 
stant and spatially isothermal gas density distribution has been 
described as this fits best the observed period distribution and 
seems a realistic assumption for the here studied case of an 
ONC-type cluster. 

In the following an isothermal gas density distribution with 
Pmax - lO^cm"^ is used for the example of an ONC-like cluster 



sound speed 


0.15 km/s 


0.49 km/s 


1.54 km/s 


maximum density 








10'*cm-^ 


37.4% 


2.8% 


0.1% 


lO'cm'' 


54.9% 


12.0% 


0.2% 


lO^cm-3 


69.2% 


28.7% 


1.2% 



Table 1. Percentage of merged binary systems for different 
sound speeds and maximum densities after 1 Myr 



In Fig. |2]the number of merged binaries Nm relative to the 
initial number of binaries N^ot, N„, - N,„/Ntot is plotted as a 
function of the initial semi-major axis of the merged binary sys- 
tems. The binaries with small semi-major axis merge to a higher 




1 10 
semi-major axis [AU] 

Fig. 2. The number of merged binaries relative to the initial num- 
ber of binaries after 0.01 Myr (dashed), 0.1 Myr (dotted) and 
1 Myr (solid) for a binary distribution embedded in an r^- gas 
density distribution with p,„aj. = 10^ cm"-' and a sound speed of 
0.49 km/s. 



degree than wider binaries. Whereas binaries with initial semi- 
major axis of 10 AU rarely merge, 40% of all binaries with initial 
semi-major axis of 0.1 AU coalesce after 1 Myr (solid line). 

The resulting semi-major axis change depends as well 
strongly on the system mass Mjys. Figure |3] shows the mean 
relative semi-major axis change Sgq = [fl(0) - a{t)]/a(Q) as a 
function of the system mass. This illustrates that the semi-major 
axis of binaries with a high system mass shrinks faster and to 
a larger degree than the semi-major axis of binaries of a lower 
system mass. The shorter coalescence time for massive systems 
(see Eq. |4|i results from the increased angular momentum and 
energy loss of high-mass binaries due to dynamical friction (Eq. 
51 in Stabler ( 2010| l). The average separation between binary 



systems with Msys > 10 M© dwindles to less than half its ini- 
tial value after 1 Myr. Additionally, the location of the binary 
system inside the cluster influences the coalescence time. For 
a maximum density of 10^ cm"-' in the cluster core, the density 
drops to 6.4 10^ cm^^ at 2.5 pc. Therefore the coalescence time 
of a binary system with Msys = lOM© drops from 0.23 Myr in 
the cluster center to 35.95 Myr at 2.5 pc. Considering this, the 
effect of orbital decay could even boost the mass segregation of 
the cluster. 

The overall effect of the gas -induced orbital decay on the 
period distribution of the binaries in the cluster is to reshape it 
by pushing binaries to tighter orbits which can even lead to the 
merging of a binary. This can be seen in Fig.|4] which shows the 
relative number of binary systems Nt - Nh{t)/NbiO) as a func- 
tion of the semi-major axis after Myr (thin solid line), after 
0.01 Myr (dashed Une), 0.1 Myr (dotted Hne) and 1 Myr (thick 
solid line) for a binary distribution embedded in a r"^ gas den- 



sity distribution with p„ 



10 cm and a sound speed of 



0.49 km/s. Because the orbital decay acts faster for tighter bina- 
ries (see. Eq.|4]i, binaries in the left part of the period semi-major 
axis distribution are altered first leading to a depopulation of bi- 
naries with semi-major axes between 0.02AU and lAU and the 
formation of a tail for semi-major axis less than 0.02AU after 
O.OlMyr. As time goes on, binaries with even larger orbits are 
affected until the orbital decay process stops when the gas is ex- 
pelled from the cluster At this point the process of orbital decay 
alone is responsible for the binary frequency in the cluster to 
drop from it's initial value of 100% to 88% through the merging 
of very tight binaries. 
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1 10 
system mass [M-,] 



Fig. 3. Mean relative semi-major axis change as a function of the 
system mass for a binary distribution embedded in an r^^ gas 
density distribution with p„,ax = 10^ cm"^ and a sound speed of 
0.49 km/s after 0.01 Myr (dashed), 0.1 Myr (dotted) and 1 Myr 
(solid). 
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Fig. 4. Chosen birth semi-major axis distribution (thin solid line) 
and effect of the gas induced orbital decay on a semi-major axis 
distributions embedded in an r^^ gas density distribution with 
Pmax -10^ cnT^ and a sound speed of 0.49 km/s after 0.01 Myr 
(dashed), 0.1 Myr (dotted) and 1 Myr (thick solid). 



4. Comparison with observations 

In the following we compare above described results with ob- 
servations. In order to do so, we converted our semi-major axis 
distributions into period distributions and re-binned our data ac- 
cording to Raghavan et al. ( 2010| l. Their observations were lim- 
ited to F6 to K3 dwarf stars which roughly corresponds to binary 
systems with primary masses Mji e [0.5 Mq, 1.5 Mq]. 

Figure 5^) shows the resulting initial period distributions as 
thin solid lines and the final period distributions as thick solid 
lines for the process of orbital decay. In accordance to the results 
for the semi-major axis distribution the orbital decay (Sec.|3]l re- 
duces the number of short period binaries. In above model the 
influence of the gas on the binary system has been studied, but 
the dynamical interactions between the binaries with other stel- 
lar components in the cluster were neglected. It is known that 
dynamical interactions destroy long period binaries. The ques- 
tion arises whether there is any region in the period distribution 
which is significantly affected by both processes? 

To investiage this, we want to show the pure dynamical im- 
pact of the cluster evolution on the binary population in absence 




lO-" 10" 
period [d] 



0.14 



b) 




10' 10" 10" 
period [d] 



10" 



Fig. 5. a) Chosen birth period distribution (thin soUd Une) and 
effect of the gas induced orbital decay on a period distribu- 
tions embedded in an r^^ gas density distribution with p,„a^ = 
10^ cm"^ and a sound speed of 0.49 km/s after 1 Myr (thick 
solid line) compared with the Gaussian fit (dotted line) to the 
observations by [Raghavan et alT] ( |2010) . b) Same as in a), with 
results from the nbody6 simulations. 



of the "gas-effect" described above. So all clusters have been 
set up without any gas and are initially in virial equilibrium. 
Only the first 3 Myr are investigated in this study. As only the 
most massive stars evolve significantly in this time, stellar evolu- 
tion was not considered in the simulations. To improve statistics, 
50 realisations of these clusters have been simulated and their 
results have been averaged. For more details on the numerical 
method see |Kaczmarek et aLl ( |201 l| l. 

To provide comparable initial conditions as in Sec. [3] all 
simulations have been set up as described in Sec. [2] The for- 
mal initial binary frequency has been chosen to be 100% i.e. 
each star is generated as part of a binary system. However, the 
resulting initial binary frequency is always lower as a conse- 
quence of the setup procedure. All particles in the simulations 
are created as binaries initially with a semi-major axis sampled 
from the log-uniform distribution. Afterwards these binaries are 
placed randomly in the cluster independent of the semi-major 
axis. Therefore, binaries with separations longer than the local 
mean separation will often have a closer partner already, so that 
thereby are fewer large-separation binaries than expected form 
the set up. 

Figure |5})) shows the period distribution of the simulated 
cluster initially (thin solid line) and at the age of 3 Myr (solid 
line) with the Gaussian fit (dotted line) to the observations. It 
can be seen that dynamical evolution of the binary population 
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in a star cluster mainly destroys wide binaries, but does not af- 
fect tight binaries. These wide binaries are destroyed in three- 
and four-body encounters where a perturber transfers some of 
it's kinetic energy to the binary system so that after the interac- 
tion, the binding energy of the binary becomes positive (IHeggie 
1975\ . However, this process is only effective if the binary is not 
too strongly bound which implies that predominantly wide bina- 
ries are affected. Nevertheless, after 3 Myr still a small amount 
of wide binaries exists. Our results give no definite answer to the 
origin of these wide binaries. They are either remainders of the 
initial binary population or formed by dynamical processes. The 
reason is that also we perform a large number of simulations the 
results for the wide binaries still suffer from low number statis- 
tics. 

Similar results were obtained by e.g. 'KroupaJ ( 1995a i; 
Kroupa et al.| (2001 ) using a initial period distributions with a 
higher number of initial wide binaries than used here (dashed 
line in Fig.[T]l- It should be noted that, using a log-normal distri- 
bution with a initial lower number of wide binaries (thick grey 
line in Fig. [T) Par ker et aL] ( [2009) obtained a deficiency of wide 
binaries in comparison to the field. If the initial binary distribu- 
tion would be log-normal, this would require wide binaries to be 
formed by other processes see for example [Kouwenhoven et aT] 
( |2010| l; [Moeckel & Batel ( |2010| l. However, apart from the wide 
binaries, the overall result does not seem to be very sensitive to 
the initial distribution. In our work a log-uniform period distri- 
bution as for the investigation of the orbital decay is applied. 
Otherwise the treatment is very similar to previous work. 

The destruction of the wide binary systems is accompanied 
by a reduction of the binary frequency by 13% + 2% after 3 Myr 
where the error is given by the standard deviation among the 50 
realisations. However, this represents only a lower limit for the 
destruction of binaries by dynamical interactions. As has been 
mentioned before, although all stars are intended to be part of a 
binary initially, a certain fraction of all binaries - dominantly the 
wide - is 'disrupted' during the setup process. Treating these as 
being destroyed by the cluster dynamics results in an upper limit 
of 26% ± 2%. In summary, the eff'ect of the dynamical destruc- 
tion of wide binaries is a reduction of a; 13% - 26% of the binary 
frequency during the first 3Myr. Comparing this to the 12% loss 
of binaries by merger caused by orbital decay, this means that 
both processes are of equal importance, at least for the consid- 
ered case. 

In the following we specify the maximum effective period for 
the orbital decay f orb and the minimum effective period for the 
dynamical destruction f dyn at which at least 3% of the binaries 
are affected by the corresponding processes (see Fig.|6]l. The re- 
sulting values are f orb - 5.5x10^ days and Pdyn - 1.1x10^ days, 
it follows Porb < ^'dyn- So there is no period range where both 
processes simultaneously play a role. Separating the two process 
as done here, does not take into account binaries that exchange 
partners and become more strongly bound. However, testing this 
for the here considered first 3 Myr of the cluster development, 
we found that only in <K 1% of all cases binary harding leads to 
a transgression into the regime where orbital decay takes place. 
This means that both processes can be treated separately (see 
Fig. |6| in a additive way. 

In order to determine the likelihood that our simulated period 
distributions correspond to the observed field popu lation we per 
formed ;^f^-tests as described by 



Press et al. 



(2007 1: 



where is the number of bins, Rj and 5, are the two distribu- 
tions to be tested. This is only applicable if J] - YjS i- To war- 

rant this we scaled our probability distributions to the total num- 
ber of binaries observed by Raghavan et al.J( |2()10 i. Afterwards 
the -probability function Q(x^\v) is used to calculate the prob- 
ability that both distributions have been sampled from the same 
underlying distribution. 



Table |2| shows the results of;\f^-tests for the Raghavan et al. 
(2010) period distribution and the simulated one. Obviously, the 
probability that the initial period distribution and the observed 
period distribution of the field have been sampled from the same 
underlying distribution is negligible (probabilities of 3.5 x 10"^ 
and 4.0 x lO^^*). Similarly, considering only one of the two pro- 
cesses "orbital decay" or "dynamical interaction" alone results 
in (higher, but still) low probabilities. Restricting the ;if^-test to 
the proper period ranges (10 ' - lO^d for the orbital decay and 
10^ - 10'°d for the cluster dynamics) yields much better re- 
sults - both the orbital decay and the cluster dynamics result 
in Q(x^\v) values of 67% and 92%, respectively, indicating al- 
ready the really good agreement of our simulated period distri- 
butions and the observations by Raghavan et al. ( 2010[ ). A su- 
perposition of the distribution of periods P < 10 d depleted by 
orbital decay and the distribution of periods P > lO^d resulting 
from the dynamical destruction is shown as thick solid line in 
Fig.[7] Performing the;^f^-test on this complete distribution over 
the entire period range results in a e(rV)-value of 94.1%. This 
clearly demonstrates, that together these two processes naturally 
reshape a log-uniform distribution to a log-normal distribution 
as observed in the field today. 





orbital decay 


N-body dynamics 






complete period range 








x' QixM 


initial 


71.5 


2.3 X 10-" 


23.6 0.009 


final 


36.4 


7.2 X 10-5 


22.4 0.013 










adopted period ranges 




x' 




x^ QixM 


initial 


45.7 


1.1 x 10->* 


6.6 0.171 


final 


3.2 


0.67 


0.9 0.92 



Table 2.;if^-test results calculated for the initial and final theoret- 



ical period distributions and those observed by Raghavan et al. 
( |2010| l. The first four column block shows the results of tests 
against the complete Ragha van et al.| ( |2010 1 period distribution 
ranging from xlO-' - xlO'" d. The second four column block 



shows tests against the left part of the Raghavan et al. 
period distribution (xlO-^ - xlO^d) for the orbital decay distri- 



( |20T0l l 



butions and the right part of the jRaghavan et aL (20101 period 
distribution (xlO-^ - xlO'^d) 



2 _ ■<r^ {R, - S if 
^ -Zj R. + Si 



(9) 



5. Conclusions 

In this paper we demonstrated for the first time that gas-induced 
orbital decay of binaries in the embedded cluster phase signifi- 
cantly changes the properties of short -period binary distribution. 
It turns out, this process is of equal importance for the develop- 
ment of the initial binary population as the well-studied process 
of dynamical evolution, at least for ONC-like clusters. 

Gas-induced orbital decay not only changes the binary prop- 
erties, but can even cause mergers of binaries, creating a more 
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1.5 
time [Myr] 

Fig. 6. Evolution of P^y„ (upper line) and f orb in the first 3 Myr 
of cluster evolution. 




10-^ ICH 10" 
period [d] 



Fig. 7. Comparison of the period distributions resulting from 
the orbital decay of embedded binaries and the dynamical de- 



struction (thick solid line) with the observations by Raghavan 
et al.|([2010) (d ashed line). Additionally, the log-normal fit by 
Raghavan et al. (2010 1 is shown as short-dashed line. 



massive single star The likelihood of such mergers depends on 
the density profile of the gas and on the separation and the mass 
of the binary system. For example, while binaries with semi- 
major axis of lOAU rarely merge after 1 Myr in a cluster with 
a maximum gas density of 10^ cm"-', 40% of all binaries with a 
semi-major axis of 0.1 AU do so. In general the separation be- 
tween binary systems with Msys > lOM© shrinks on average to 
less than half its initial value after 1 Myr. 

Orbital decay changes only the period distribution of short 
period binaries. However, long-period binaries, are affected by 
dynamical destruction caused by the interaction of the stars in a 
dense cluster environment. Our simulations show that interme- 
diate period binaries are nearly unaffected by either of these two 
processes (see Fig. [8]). 

For an ONC-like cluster we found that for these G-type pri- 
mary stars that the orbital decay due to the interaction of the 
binaries with the ambient gas basically only affects binaries with 
periods smaller than 5 x lO'^d corresponding to separations closer 
than X 36AU. By contrast, dynamical interactions destroy only 
binaries with periods larger than lO^d (i.e. x 53AU). This means 
that there is no region in the period distribution that is influenced 
by both effects. Therefore it is possible to investigate the two 
processes separately and combine the results. 



orbital decay 



V 

1 



initial 

cluster dynamics 



final 



10" 



10' 



10-' 



10' 



10' 



10» 

period [days] 



Fig. 8. Schematic illustration of the two-stage process forming 
the log-normal field-binary period distribution observed today 
(grey line) from an initial log-uniform distribution (solid line). 



Perhaps the most striking result of this investigation is, 
that these two processes transform a primordial log-uniform 
into a log-normal period distribution, on a realistic time scale. 
The distribution resembles, in all its properties, remarkably 
that observed in the field, without need of further assumptions. 
Performing a;(f^-test of the resulting distribution and recent ob- 
servations of the field binary population ( Raghavan et al.|2010| l 
yields a probability of 94.1% that both distributions originate 
from the same origin. 

The emerging picture can be summarized the following way 
(illustrated in Fig. [8]): Binaries bom in dense, embedded clusters 
with initially log-uniformly period distributions are processed 
in two ways. As long as the cluster is embedded in its natal gas 
(about IMyr), the orbital decay of the embedded binaries depop- 
ulates the left hand side of the period distribution. The dynamical 
evolution of the cluster destroys wide binaries, depopulating the 
right hand side of the period distribution. The combined effect 
of these equally important processes is that the final period dis- 
tribution of the binary population in the star cluster has become 
log-normal although it initially has been log-uniform. 

This re-shaping of the period distribution leads to a reduc- 
tion of the binary frequency by 29%. The gas-induced orbital 
decay shrinks the separation between two stars such that merg- 
ers in 14% of all cases are possible. The dynamical interactions 
destroy 15% of G-type binaries by three-body encounters. This 
latter fraction only represents a lower limit to the reduction of the 
binary frequency due to out specific setup procedure and can be 
as high as 28%. Consequently the overall reduction of the binary 
frequency could be up to 42%. Assuming that most stars form in 
an ONC-like cluster and adopting a field binary frequency of G- 
type primaries is 58% (Ragh avan et al.|2()T0 ) our results strongly 
favour a much higher initial binary frequency of 87% - 100% 
(assuming 44% one gets 73 - 86%) for solar-type stars. 

Mass segregation and high densities in the cluster centre 
favour the merging of massive stars. From the merging of pre- 
dominantly massive stars, one would expect to see a difference 
in the IMF of single and binary stars, with the single star distri- 
bution having an excess of massive stars. This has not been ob- 
served. However, the situation is more complex. Massive stars 
are as well the most likely to capture a new partner For exam- 
ple mergers which form a stars with M > 20M|T[ , can happen 
on a timescales of lQ*yr (jPfalzner & Olczak 2007 1. The merger 
product and its companion would become part of the binary pop- 
ulation again. So the IMF of the single and binary population 
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would be changed by merger processes. Weather they would be 
effected to the same degree would require further studies, but 
given the poor statistics at the high mass end, they would be in- 
distinguishable. 

Clearly, the field distribution is a mixture of stars originating 
from stars that experienced star formation in isolation, sparce 
(Taurus-like), dense (ONC-like) and very dense (Arches-like) 
clusters. On the one hand, associations like the Taurus clusters 
probably never had a gas density above 1 0"^ cm- 3 and the cur - 
rent stellar density is below 10 stars pc^^ ( Luhman et al.|2069 1. 
This means in such systems the binary population will neither be 
affected by orbital decay nor orbital destruction in a significant 
way and is therefore nearly unaltered from its primordial state. 
On the other hand, orbital decay becomes important for high gas 
densities and dynamical destruction for high stellar densities 

On the basis of actual observations, one can only speculate 
how many stars are bom in these different density regions. While 
Bressert et al.| ( [20^1 0) stated, that only a minor part of all stars 
in the solar neighbourhood form in high density regions, this is 
unclear for the total of the galaxy. |Dukes & KrumhoIz| ( |201 l| l 
concluded, that 1/2 - 2/3 of all stars are in clusters with more 
than 1000 stars. As such massive clusters initially had a much 
higher stellar densities ( |Pfalzner|2009 ), this indicates that envi- 
ronmental effects are important for the latter The remarkable re- 
semblance of the distribution, obtained after both processes have 
taken place, to the log-normal period distribution of the field bi- 
nary population is a strong indication, that ONC-like clusters 
might be a dominant contributor to the field distribution. 

Additional to the here studied ONC-type clusters, further in- 
vestigations should include a variety of initial conditions: differ- 
ent stellar and gas density distributions (Kroupa,, 1995b, ,Parker, 
et al. 201 l|i, a range of cluster densities (Olczak et al. 2010 



Marks et al.||201l l and different virial states of the cluster 



(Allison 2009). Currently we test the limitations of the model 
of orbital decay by numerical simulations (Korntreff & Pfalzner, 
in prep). 

To verify our initial conditions observationally, detailed stud- 
ies of binary populations in very young (< 100000 yrs) embed- 
ded star clusters, before the onset of the here investigated pro- 
cesses, would be necessary. However conclusions about the bi- 
nary distribution in such young systems would be hindered by 
low-number statistics. Usually only a few tens of stars are ob- 
servable, although the true membership might be considerably 
higher due to the high extinction in these clusters. Combining 
data from different such clusters would be an alternative. An 
other approach could be to observe different regions within a 
single OB associations. The processes are predominant at work 
in the association centre and much less so in the outskirts. 
Therefore, it could be expected that the period distribution in the 
association centre differs considerably from the outskirts - close 
to the log-normal in the cluster center and log-uniform in the 
outskirts. However, this picture neglects the mixing within the 
cluster, so that further theoretical work is needed to investigate 
this point. 
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